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Abstract. We consider quantum Heisenberg ferro- and antiferromagnets on the square lattice with exchange 
anisotropy of easy-plane or easy-axis type. The thermodynamics and the critical behaviour of the models 
are studied by the pure-quantum self-consistent harmonic approximation, in order to evaluate the spin and 
anisotropy dependence of the critical temperatures. Results for thermodynamic quantities are reported and 
comparison with experimental and numerical simulation data is made. The obtained results allow us to 
draw a general picture of the subject and, in particular, to estimate the value of the critical temperature 
for any model belonging to the considered class. 

PACS. 75.10.Jm Quantized spin models - 75.10.-b General theory and models of magnetic ordering - 
75.40.Cx Static properties 



1 Introduction 

In this paper we examine the class of magnetic models 
described by the Hamiltonian 

H = -jJ2 (S?S? + Sf S*j) + ASfSf] (1) 

<ij> 

where i = (ii 5*2)3 j = (jit 32) & r e sites on a square lattice, 
and the sum runs over the pairs of nearest neighbours; 
the spin degrees of freedom are described by the quan- 
tum operators Si = {Sf, Sf , a"), obeying the angular mo- 
mentum commutation relations [Sf, <Sj- ] = e a ^ 1 8\- } S7 with 

|S| 2 = 5(5+1). 

For positive A and /Lt, the sign of the exchange inte- 
gral J sets the ferro- or antiferromagnetic character of the 
model, J > or J < respectively; for the sake of a more 
compact notation, however, we will hereafter take a posi- 
tive J and properly change the sign of fi and A to obtain 
the required antiferromagnetic Hamiltonian. 

At a classical level, i.e. when 5 — > 00, ferro- and anti- 
ferromagnets have the same static properties and can be 
simultaneously studied as far as their thermodynamics is 
concerned. On the other hand, when quantum fluctuations 
are taken into account, the two types of models display 
substantially different features. It is worthwhile noticing 
that, despite displaying the most interesting quantum fea- 



tures, antiferromagnetic models are less affected by quan- 
tum fluctuations than ferromagnetic ones, as quantum 
renormalizations are sensibly stronger in the latter, what- 
ever the approach used to evaluate them. 

In this paper, the anisotropy parameters /i and A are 
set to range in the interval [—1,1], and the following sub- 
classes of models are hence considered: 

(i) Isotropic: /1 = A, |A| = 1, ferro- and antiferromag- 
netic (A > and A < 0, respectively). 

(ii) Easy-plane: /1 = 1, A 6 (— 1, 1), ferro- and antifer- 
romagnetic (A > and A < 0, respectively). 

(hi) Easy-axis: fi G (—1,1), A = sign(/i), ferro- and 
antiferromagnetic (/1 > and \i < 0, respectively). 

Models belonging to these sub-classes have the com- 
mon feature of displaying, as the temperature decreases, a 
phase-transition towards a more ordered phase; this tran- 
sition generically occurs at a critical temperature which is 
a function of both the spin value and the anisotropy pa- 
rameters, and will be hereafter indicated by tc^is, 5), with 
v = fx,X depending on the specific class considered, and 
the reduced temperature t defined by t = T/JS 2 with 
5 = 5 + 1/2. 

(i) The isotropic model, both ferro- and antiferromag- 
netic, has the full rotational 0(3) symmetry and it is hence 
in a paramagnetic phase down to t = 0, as required by the 
Mermin- Wagner theorem Its ground state, however 
and at variance with the one-dimensional case, is ordered 
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for all spin values also in the antiferromagnetic case Q, 
and the model can hence be thought to have a phase tran- 
sition at a critical temperature t c (H,S) = V5*. 

(ii) In the easy plane case, the Mermin- Wagner the- 
orem still holds and no finite temperature transition to- 
wards a phase with a finite order parameter, may occur. 
However, a Berezinskii-Kosterlitz-Thouless (BKT) transi- 
tion HQ, related with the existence of vortex-like topo- 
logical excitations, is known to characterize the class of 
the easy-plane models and may occur at a critical tem- 
perature £bkt(A, S) > 0. 

The reference system for the easy-plane class is the pla- 
nar, or XY, model, defined by Eq.(|) with A = 0, fi = 1. 
Its BKT critical temperature has been seen to be finite for 
all spin values: above £bkt the system is disordered, with 
exponentially decaying correlation functions. In the whole 
region < t < £bkt the system is in a critical phase with 
vanishing magnetization: its correlation functions decay 
according to a power law, testifying to the existence of 
gwasz-long-range order; at t = the magnetization gets a 
finite value and the system is ordered. The XY model may 
be thought to display two phase-transitions: a BKT one 
at a finite temperature, followed by one towards a fully 
ordered phase at t = 0. 

(iii) In the easy-axis case, the Mermin- Wagner theorem 
does not hold, and the easy-axis models may undergo a 
transition of the Ising (I) type towards a phase with long- 
range order, at a critical temperature t\(n, S) > 0. 

The reference system for this class of models will be 
hereafter called the Z model, described by Eq.(Q) with 
fi = 0, A = 1; such model has a transition temperature 
ii(0, S) > for all spin values. Above t\ the system is in 
a paramagnetic phase, with exponentially decaying cor- 
relation functions; below t\ the system is ordered with a 
finite magnetization along the easy axis. The Z model 
with S = 1/2 is the Ising model, representing an impor- 
tant source of information for the whole class of the easy- 
axis models, not only because of the existing exact results 
by Onsager ||, but also because of its fundamental role 
in the renormalization-group approach and in conformal 
field theory. 

The information on the reference models, as well as 
the general considerations based on the symmetries of the 
Hamiltonian, are surely valuable, but not sufficient to fully 
appreciate the richness of features contained in Eq. ([!]), 
and to properly understand the thermodynamic and criti- 
cal behaviour of the many existing real compounds whose 
magnetic behaviour is described by such Hamiltonian. To 
accomplish this goal one has to study how the specific 
values of A, fj,, and S may determine the behaviour of the 
system. 

At a qualitative level, we know that quantum fluctu- 
ations, whose strength is measured by the quantum cou- 
pling 1/5, introduce a disordering agent in the classical 
(S —>■ oo) picture: as a consequence, we expect the possible 
critical temperatures to decrease with decreasing S. The 
same effect is obtained by increasing |A| (easy-plane case), 
or (easy- axis case): this in fact means to weaken the 
anisotropic character of the system, thus allowing larger 



spin fluctuations out of the easy plane or the easy axis, 
and inducing a lowering of the critical temperature. To 
this respect it is to be noticed that most real compounds 
have very weak anisotropies Q, i.e. values of A or /i very 
close to unity: this is the region where quantitative predic- 
tions are more difficult to obtain, as the model's character 
is not well defined and fluctuations are large. 

Aim of this paper is to provide a quantitative descrip- 
tion of the model (|l|) when S and A, or [i, are varied. 
We will mainly concentrate on the values of the criti- 
cal temperatures, but also give expressions and results 
for other thermodynamic quantities. The method used 
is the pure-quantum self-consistent harmonic approxima- 
tion (PQSCHA), extensively described in Ref. (7| and al- 
ready applied to many magnetic systems (see for instance 
Rcfs. H||). Its implementation is markedly different de- 
pending on the specific class of models considered, as ex- 
plained in Section |^, and the easy-plane and easy-axis 
cases will be hence considered separately (in Section || 
and |, respectively) . In Section || we will comment on the 
resulting global picture and compare our results with all 
the available experimental and numerical simulation data. 
Conclusions about the critical behaviour of the models will 
be finally drawn in Section o. 



2 Method and spin-boson transformations 

The PQSCHA is a semiclassical method based on the 
path-integral formulation of quantum statistical mechan- 
ics: the main peculiarity of the method is that of defin- 
ing an analytical separation between classical and pure- 
quantum contributions to the thermodynamics of the sys- 
tem jH^, thus allowing the possible exact treatment of 
the former, while requiring a self-consistent harmonic ap- 
proximation of the latter. This means that the only ap- 
proximated contribution is the non-linear pure-quantum 
one, which is in fact considered at the one-loop level. The 
value of such result is evident if one considers that most of 
the peculiar features of magnetic systems, including their 
critical behaviour, are determined by long-wavelength ex- 
citations, whose character is essentially classical; the pos- 
sibility to fully take into account their non-linearity at 
the classical level, together with the fact that harmonic 
pure-quantum fluctuations are exactly considered, justi- 
fies the success of the PQSCHA in predicting most ther- 
modynamic properties of magnetic systems in a large tem- 
perature range. 

Without going into the details of the derivation of the 
PQSCHA, we report here the essential formulas, and dis- 
cuss the scheme of implementation to magnetic systems, 
where differences between the sub-classes of models be- 
come essential. 

Consider a quantum spin system on a lattice described 
by the Hamiltonian H(S) with S = (Si,...Sjv) and N 
number of lattice sites; if O is the quantum operator rela- 
tive to a physical observable of the system, the PQSCHA 
expression of its quantum statistical average has the clas- 
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sical-like form 



:°> = -5- / d N sO ee (s)e 

ArfF 



the symmetry properties of the original Hamiltonian and 
of its ground state. 

-/3W e ff(s) (2) I n markedly easy-plane cases (|A| <C 1) we use, for each 

spin operator, the Villain transformation (VT) JL2| 



where (3 = t 



--off 



Id' 



se 



-/3H 



eir 



and s = (si, ...sn) 



with Si classical unit vectors. The effective Hamiltonian 
H c ff(s;f, S), is a classical-like function, depending on the 
parameters of the original Hamiltonian, and determined 
by the PQSCHA renormalization procedure; a similar pro- 
cedure, when applied to the quantum operator O, leads to 
the function O e g(s; t, S) 

The phase-space integral in Eq. (|^) may be evaluated 
by any classical technique, such as the classical Monte 
Carlo (MC) simulation. Eq. (J2J) , on the other hand, con- 
tains quantum renormalizations embodied in the temper- 
ature and spin dependence of Ti c s and o ff ; such depen- 
dence causes the effective classical model to change for 
each temperature point, so that ad hoc simulations must 
be performed, and only in very peculiar cases one can di- 
rectly use existing MC data (see for instance the case of 
the isotropic model in Ref. ||). 

The PQSCHA naturally applies to bosonic systems, 
whose Hamiltonian is written in terms of conjugate oper- 
ators q = (qi, ...Qn), P = {pi, -Pn) such that [q m ,p n ] = 

i5 mn ] the method, however, does not require Ti.(p, q) to be 
standard, i.e. with separate quadratic kinetic p-dependent 
and potential g-dependent terms, and its application may 
be extended also to magnetic systems, according to the fol- 



lowing scheme |10 : The spin Hamiltonian 7i(S) is mapped 
to 7i(p, q) by a suitable spin-boson transformation; once 
the corresponding Weyl symbol H(p, <?), withp = (pi, ...pn) 
and q = (qi,...qjsr) classical phase-space variables, has 
been determined JTIJ , the PQSCHA renormalizations may 
be evaluated and H c s{p,q) and O c s{p,q) follow. Finally 
7i ff(s) and G ff(s) are constructed by the inverse of the 
classical analogue of the spin-boson transformation used 
at the beginning. 

In order to succesfully carry out such renormalization 
scheme, the Weyl symbol of the bosonic Hamiltonian must 
be a well-behaved function in the whole phase space. Spin- 
boson transformations, on the other hand, can introduce 
singularities as a consequence of the topological impossi- 
bility of a global mapping of a spherical phase space into 
a flat one. The choice of the transformation must then be 
such that the singularities occur for configurations which 
are not thcrmodynamically relevant, and whose contri- 
bution may be hence approximated. Most of the meth- 
ods for studying magnetic systems do in fact share this 
problem with the PQSCHA; what makes the difference 
is that by using the PQSCHA one separates the classical 
from the pure-quantum contribution to the thermal fluc- 
tuations, and the approximation only regards the latter, 
being the former exactly taken into account when the ef- 
fective Hamiltonian is cast in the form of a classical spin 
Hamiltonian. 

The spin-boson transformation which constitutes the 
first step of the magnetic PQSCHA is chosen according to 



S~ 



0$+)+, S z =p; (3) 



this transformation keeps the 0(2) symmetry in the easy 
plane, meanwhile allowing to deal with the square root in 
terms of a physically sound small-p approximation. 

In the easy- axis case, it makes no sense to use Eqs.(||), 
as the expectation value of the z-component of each spin is 
now substantially different from zero. On the other hand, 
the Holstein-Primakoff transformation (HPT) |Q 



V2S+ = \jS + S-z 2 (q + ip) , 

= (£+)*, S z =S-z 2 , (4) 

and the Dyson-Maleev one (DMT) JTJ] 



V2S + =q + l P , S- = (q - ip)(S + S ~ z 2 ) , 



(5) 



with z 2 = (q 2 +p 2 )/2, do both suggest the z-component 
to be the privileged one for alignment, thus fitting to the 
\fi\ < 1 case. 

The isotropic, or nearly isotropic, case is the most dif- 
ficult to be treated, as no part of the spherical phase-space 
may be generically chosen to be described better than an- 
other; nevertheless, when the ground state is at least char- 
acterized by long-range order with an alignment axis, the 
HP and DM transformations are both seen to work well 
in the PQSCHA scheme. The reason is that the boson 
picture only serves to derive the pure-quantum renormal- 
izations, while the symmetry is thereafter restored in 7i c ff 
so that it is sufficient to assume a local alignment axis. 



3 The easy-plane case 

The Hamiltonian Eq. (0) in the easy-plane case takes the 
form 



n 



-j 



<ij> 



as? sy 



(G) 



with A € (— 1, 1). The antiferromagnetic sector is mapped 
in the A € (—1,0] region of Eq. (^) by the canonical trans- 
formation (S?,S?,S?) -> ((-l) l Sf, (-lySf^f) where 

(— l) 1 = (— 1)' 1+ * 2 ; ferro- and antiferromagnets may be 
thus simultaneously studied. 

The effective Hamiltonian relative to Eq. ([|) is 



Hes = —JS 2 j, 



eff / 
<ij> 



S\S V A +A eff sf Sf 



(7) 
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where the renormalization parameters j e s(t,X, S), 
X e s(t, A, S), and Q(t,X,S) have different expressions de- 
pending on the A-region considered. 

In the A< 1 case, the VT Eqs. (||) may be safely used, 
and give pE| 



Jeff 



1 



= A(1- 

g v = js 2 n 



D 



) 2 e - C ll/ 2 



A 



tr-im(i-\D,)-^ 



(8) 
(9) 
(10) 



where the renormalization coefficients D± and X>|| repre- 
sent the effects of in-plane and out-of-plane pure-quantum 
fluctuations, respectively, and have the form 



1 



2SN ^ a k 

J^V(l- 7k )^£ k 
2SN k 6k 



with 



4 e - c "/ 2 (l- A off7k ) , 
4(1 



)2 e -^|/2 ( 



l-7k) 



(11) 
(12) 

(13) 
(14) 



£ k = (coth/ k - l// k ), r = N- 1 E k ln(sinh/ k // k ), / k = 

a k b k /(2St), 7k = £ d cos(k ' d )/ 4 ' d = (±1.±1)> and k 
is a wave vector in the first Brillouin zone. The uniform 
coefficient <9 V is 



-fn/2 (2D ± + (1 - D_l)D| 



Close to the isotropic limit, A < 1, the HPT must be 
used: because of the residual easy-plane character of the 
model, the privileged axis, i.e. the quantization one, must 
be chosen in the easy plane; in particular, we choose the y 
direction as the local alignment one. A quite complicated 
PQSCHA implementation finally leads to the following 
expressions jL6| : 



Jeff ^ 



G HP = JS' Z N 



d 



npp 

(l-X) + —(l-X) 



tr- 



(15) 
(16) 

, (17) 



where the temperature and spin dependent renormaliza- 
tion coefficients are defined as follows: 

2 = 1 - \v + \d , d = - X cS D pp ) , (18) 



V 



— Ed-7k)^ k , 

k 



and 
D pp 



2SN ^ W 



D qq 



2SN ^ a k 



(19) 



7k £k , 



1.0 



'eff 
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Fig. 1. Renormalization parameter j c ff vs t for the easy-plane 
model with A = (full curve), 0.5 (dashed), 0.7 (long-dashed), 
0.9 (dash-dotted), and S = 1/2, 1, 5/2. The dotted line is j cff 
for S = 1, X = 0.5, as obtained by the HPT. 



with dimensionless frequencies 



40 2 (1-A cff7k ) , 
b 2 k = 40 2 (l- 7k ) . 



(20) 
(21) 



The uniform coefficient (9 HP is 



hp = Dqq ^ _ 3D q _ D pj + xd pp( 3 _ 3D p _ D g) + 

+ 2{D qq ) 2 + 2{D PP ) 2 - (D q + D p ) 2 + 2{D q + D p ) , 



with 



D p = 



2SN 



&k 



D q = 



2SN ^ a k 



The renormalization of the exchange energy by the fac- 
tor j c g < 1, and the weakening of the easy-plane anisotropy 
(|A| < |A c ff|), are seen to be, both for strong and weak 
anisotropy, the result of the cooperative effect of in-plane 
and out-of plane pure-quantum fluctuations. It is worth- 
while noticing that, in the isotropic limit |A| = 1, the 
effective Hamiltonian (0), determined by the HPT, coin- 
cides with the one found directly for the isotropic model 
by the DMT. 

In Fig. [TJ we show the renormalization parameters j e g 
vs. t, for different values of S and A; the two most (least) 
anisotropic cases are treated by the VT (HPT). For 5 = 1 
and A = 0.5, j e g as obtained by the HPT is also shown 
(dotted line): the difference with the VT result (dashed 
line) is remarkable, testifying that both transformations 
cannot give very accurate results in such intermediate re- 
gion. However, as we will see in Section ||, this does not 
dramatically affect the critical temperature estimation. 

The difference (A c ff — A) is shown in Fig. || as a function 
of A, for different values of S and t. For each couple (t, S) 
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0.12 



0.08 



0.04 - 




0.00 



Fig. 2. Difference (A c ff — A) vs X for the easy-plane model with 
S = 1/2 (thin) and 1 (thick), at t = 0.5 (full) and 1 (dashed). 
For each t and S, both the HP (curves passing through (±1,0)) 
and the V ( curves passing through the origin) results are 
shown. 



the curves passing through the origin are obtained by the 
VT, while those passing through the point (—1,0) or (1,0) 
are obtained by the HPT in the antiferro- or ferromagnetic 
case, respectively; the two types of results do not smoothly 
connect to each other, but nevertheless cover the whole 
[—1,1] interval and show a behaviour of the same type 
of that found in the easy-axis case (see Fig.||), where the 
problem of using two different spin boson transformations 
is not present. As expected, the HP curves display the 
correct |A| — * 1 behaviour (i.e. |A c ff| HP — > 1 
seen to be proportional to A, when A — » 0. 

As mentioned in the Introduction, the easy-plane mod- 
els are characterized by a critical phase, with vanishing 
magnetization and diverging correlation length, extending 
in the temperature region < t < teKT! thermodynamic 
quantities that may be significantly studied in such phase 
are, for instance, the free energy and the specific heat. To 
this respect, notice that the uniform term Q, Eqs. (|To| ) and 
(|l7j), not entering the statistical averages' calculation, is 
essential in the evaluation of the partition function Z. 

In the paramagnetic phase t > £bkt, one of the fun- 
damental thermodynamic quantities to be investigated, is 
the in-plane correlation length defined by the asymptotic 
behaviour of the correlation functions 



while A^j is 



G(r) 



s 2 



s?s?_ 



(22) 



with r = (r\, r 2 ) any vector on the square lattice, r = |r| 
and r) r = 1 or (— l) r i+ r 2 j n the ferro- or antifcrromagnetic 
case. 

By using Eq. @) the following PQSCHA expression is 
found 



G(r) = A r (s 



i+r 



'off 



(23) 



where ( • ) c ff = Z^ 1 J d N s( ■ )e~° Half is the classical-like 
statistical average defined by the effective Hamiltonian. 

When the VT is used, i.e. in the strongly easy-plane 
case, Eq. (E3h holds for all values of r, with 



A 



VT 



(1 



1 



V ± ) 2 e 



-X>||/2 e -U r 



(24) 



and 



«k 



cos(k • r)£k 



In the nearly isotropic region, the use of the HPT leads 
to a much more complicated general expression for G(r), 
that nonetheless goes to expression (|23|) with coefficient 



A 



HP 



D p ) 



3D p ) 



(25) 

for large r . From Eqs. ( |24| ) and (|25j), the renormalization 
coefficient A r is seen to converge to a finite value when 
r — > oo; this property, as we will see in Section ^, is of 
fundamental importance when the critical behaviour of 
the model is under investigation. In particular it means 
that the correlation length of the quantum model is related 
with that of its classical counterpart by the relation 



(26) 



with £ e ff = tj j Q g. This relation allows us to evaluate the 
correlation length of a quantum model with spin S, at a 
given temperature t, by performing a classical MC sim- 
ulation on a model with anisotropy A e ff, at the effective 
temperature i e g; a particularly simple case is that of the 
XY model A = A c ff = 0, whose correlation length may 
be directly determined in terms of the classical one, by 
a simple temperature scaling. Fig. |^ shows the correla- 
tion length versus t, for S — 1 and A = 0, —0.5, —0.999, 
and -1. The classical MC data for A = and A = 0.999, 
Rcfs. UJI and m& respectively, are also shown. 



4 The easy-axis case 

The Hamiltonian Eq. ([[]) in the easy-axis case takes the 
form 



(27) 



<u> 



with /i 6 (—1, 1) and A = sign(/i). The antiferromagnetic 
sector is mapped in the fi £ (—1,0] region of Eq. ( |27| ) by 
the canonical transformation 

(S?,S?,S?) -> h-iySf, A. Because of the 

easy-axis character of the model, the HP and DM trans- 
formations may be safely used, and actually give the same 
results; the latter, however, may be preferred because of 
its simpler structure. 
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Fig. 3. Correlation length £ vst for the easy-plane model with 
3 = 1 and A = (full), -0.5 (dashed), -0.999 (dash-dotted), 
and —1 (l ong -dashed). Symbols are classical MC data for A = 
(squares |l7f) and A = 0.999 (triangles [[l6[). 



The PQSCHA effective Hamiltonian relative to Eq.@ 
is found to be [G~sl 



H cS = - JS 2 JcS [Meff (Asfs? + sV s|) + S ? S ? 

<ij> 

where j e s(t,fi, S), /i e ff(t, //, 5), and /«, 5) are 

Jeff — C|| , Moff — A* "TJ ' 

£ = 2JS 2 (l-6>j}) + JS 2 / 1 , 



with 



= 1 



©ii 



= 1 - 



V, 



/ sinh / k 



-a, 

(28) 

(29) 
(30) 

(31) 
(32) 



and the coefficients 



(33) 



are self-consistently determined by solving Eqs. ( pT| ) and 
(pq), with dimensionless frequencies 



al = 4 



.7k 



fc 2 = 4 (tf 2 - /.A^ 7 k 



(34) 
(35) 



As in the easy-plane case, quantum renormalizations wea- 
ken both the energy scale (J c g < 1) and the anisotropy 



Fig. 4. Renormalization parameter j e g vs t for the easy-axis 
model with /i = (full), 0.5 (dashed), 0.9 (dashed), and S = 
1/2, 1, 5/2. 



0.06 



0.04 - 



0.02 



0.00 




Fig. 5. Difference (fj, e s — fJ>) vs fi for the easy-axis model with 
S = 1/2 (thin) and 1 (thick), at t = 0.5 (full) and 1 (dashed). 



(I /•'off | > A*)- I n Fig. H we show the renormalization pa- 
rameter jcff vs. t for different values of S and /i. In Fig. |^, 
where the difference /i c ff — fi is shown as a function of//, we 
see that, at variance with the easy-plane case (see Fig. ||), 
the almost isotropic regime is smoothly connected with 
the strongly axial one, testifying to the suitableness of the 
DMT in treating the whole easy-axis subclass of models. 

When /i vanishes the renormalization coefficient T>± 
may be shown to keep finite, so that fi e g — > 0, and the 
Z model, with a single optical mode a 2 . = b 2 . = 40jj, Vk, 
is recovered; for S — 1/2 such model coincides with the 
Ising one. 

The easy-axis case is characterized by a finite tem- 
perature phase-transition towards a fully ordered phase; 



Alessandro Cuccoli et al.: Phase transitions in two-dimensional anisotropic quantum magnets 



7 



many thermodynamic properties may be hence studied 
both above and below the critical temperature t\, and, in 
particular, the PQSCHA expression for the magnetization 
along the easy axis m = J2 i T]i(Si}/NS is 



1.00 



m = OAs-^eff 



(36) 



The most significant quantity in the paramagnetic phase 
t > t\ is the correlation length, here defined by the asymp- 
totic behaviour of the isotropic correlation functions 



G(r) 



Hi 

s 2 



(Si • Su 



the PQSCHA expression for G(r) is 

G ( r ) = #r( S i ' Si+r>off , 

where 6% = 1- V r /2 and 



(37) 



(38) 



cos(k-r)) £k 



is a further site-dependent renormalization coefficient. For 
increasing r, as in the easy-plane case, the coefficient 8^ 
keeps finite and rapidly converges to a uniform term: the 
asymptotic behaviour of G{r) is hence determined by that 
of the classical-like correlation functions (si • Si +r ) c ff , and 
the relation 



£(t,/i,S) = r(ieS,Meff) 



(39) 



may be used. As in the easy-plane case, results for the 
quantum reference (Z) model are obtained by a simple 
temperature scaling of the classical data; we have hence 
performed MC simulations on the classical Z model, and 
produced results for the quantum model with S = 1/2, 1 
and 5/2. The corresponding magnetization and correlation- 
length curves are shown in Figs.^ and |7], respectively. As 
for the latter, notice the much sharper divergence with 
respect to the XY model (compare the dashed curve in 
Fig.0 with the full curve of Fig.||). 



5 Transition temperatures 

The PQSCHA expressions for the effective Hamiltonian 
and the correlation functions contain two essential pieces 
of information: (a) the effective Hamiltonian keeps the 
symmetry properties of the original quantum model; (b) 
the quantum correlation functions share the asymptotic 
behaviour with their classical counterpart. These two state- 
ments not only lead us to the relations ( p6| ) and (^) , but 
also allow us to assert that the critical behaviour of the 
quantum model is that of its effective classical counter- 
part, and that, for the critical temperature, the relation 



t c (v, S) = j c g(t c , v, S) tf(v eS (t c , v, S)) 



(40) 



holds. 

Due to the temperature dependence of j e s and u e g, 
Eq. (40) is quite involute; it may be nonetheless numeri- 
cally solved if an analytical expression for t c ^(v) is avail- 
able. Such expression does not actually exist, neither in 




Fig. 6. Magnetization curves for the Z model with S = 1/2 
(full), 1 (dashed), 5/2 (dash-dotted), and oo (dotted); symbols 
are our classical MC data. 



100 




Fig. 7. Correlation length for the Z model with S = 1/2, 1, 
5/2, and oo (curves and symbols as in previous figure). 



the easy-plane nor in the easy-axis case; we have hence 
collected a set of classical MC data properly distributed 
in the [0, 1) interval, in order to determine a precise inter- 
polating function t^{v). 

In the easy-plane case, classical MC data have been 
taken from Ref. |l7| and interpolated by the function 

TbktO) =a(l + 6A 2 + cA 4 -dln(l-A 2 )) _1 , (41) 
a = 0.695 , b = -0.012 , c = 0.032 , d = 0.0648 , 

where the A — > —A invariance has been imposed. 

In the easy-axis case, sources of the classical MC data 
have been Refs. |l!| and |2C(| , as well as our new data 
tf(0.5) = 0.88 ± 0.01 and tf(0.8) = 0.81 ± 0.01; the re- 
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suiting interpolating function is 

rf(fi) =a(l + bn 2 + cfx 1 - dln(l - 
a = 0.917, b = 0.068, c = 0.097 



, (42) 

, d = 0.0636 , 



where, again, the constraint tf(fi) = tf(—[i) has been 
imposed. 

Equations ( |4l| ) and (fl2|) should only be considered as 
reasonable interpolating functions of the available data; 
in particular, they do not reproduce the correct asymp- 
totic behaviour in the v — > 1 isotropic limit. Such asymp- 
totic behaviours, on the other hand, are not accessible 
by numerical techniques and are actually determined by 
Renormalization- Group approaches. However, all of the 
numerical and experimental data are in the ^-region where 
t c (v) is sensibly different from zero, as seen in Figs. § and 
i.e. where the above mentioned asymptotic behaviours 
have not set in yet. 

Once tf(v) ~ T^(v) has been determined, Eq. ( ^o|) 
may be numerically solved thus giving the value of the 
quantum critical temperature, for any value of v and 5, 
in a few seconds on a standard PC. 

In Fig. H we show the resulting curves for £ = 1/2, 
1, 5/2, and oo in the easy-plane case: for each spin, the 
whole (—1, 1) interval is covered by two different types of 
curves, corresponding to the use of j^ s ,X^ s and 
in the strongly easy-plane and nearly isotropic regions, 
respectively. The available quantum MC data ]2l] , ^2|| are 
also shown, and agree very well with our results, despite 
being in the region where, because of the strong quan- 
tum fluctuations, the use of the PQSCHA becomes more 
delicate; unfortunately no experimental data exist, to our 
knowledge, for £bkt- 

In Fig. H the transition temperature 5) is shown 
for 5 = 1/2, 1, 5/2, and oo; in this case, there exist several 
experimental data, for the real compounds K2MF4 (/i = 
-0.996, 5=1, ij xp = 0.416) MM, Rb 2 NiF 4 (ji = 
t exp 



1, 



-0.982, 5 

5 = 5/2, t° Kp = 
5 = 5/2, t cxp - 
5 = 5/2, 



,cxp 



, - 0.477) I25J, K 2 MnF 4 (fi = -0.995, 
0.553) |23j,|§7Rb 2 MnCl4 (p = -0.997, 
0.558) |P71, a nd Rb 2 MnF 4 (p = -0.994, 
= 0.575) p3|,p8[, whose agreement with 
our results is very good, despite being in the most difficult 
region to be studied, i.e. in the nearly isotropic one. 

The 5=1/2 easy-axis antiferromagnetic case deserves 
a few more comments, as a quite well distributed scan of 
the /j,E [— 1, 0) region may be performed by putting to- 
gether the exact result for the Ising model ||, tj(0, 1/2) = 
0.567, the quantum MC data by Ding |||, fi(0.5, 1/2) = 
0.5, ti(0.99, 1/2) = 0.298, and the experimental data for 
the real compound YBa 2 Cu 3 6 .i {p ~ 0.992, S = 1/2, 
^ x p = 0.310) |TJ. As for the PQSCHA, we know that its 
results are quantitatively reliable if quantum renormaliza- 
tions are not too strong: a reasonable criterion to check 
whether such condition is fullfilled or not is to require 
ieff ^ 0.5, meaning a quantum renormalizations of the en- 
ergy scale of about 50%. Such criterion, as seen in Fig. is 
not satisfied for 5 = 1/2 and t < 0.5: our ti(/i, 1/2) curve 
does in fact underestimate the most anisotropic numeri- 
cal and experimental data, despite displaying the correct 
qualitative behaviour. A significant improvement may be 



BKT 




Fig. 8. BKT critical temperature vs X for the easy-plane 
model with S = 1/2, 1, 5 /2 . Sy mbols are the quantum MC 
data for 5 = 1/2 (circles [ pi| , |22| ), and the classical MC data 
(squares ]l7|]) used to construct tb' kt (A), shown as dotted line. 




Fig. 9. Ising critical temperature vs fi for the easy-axis model 
with S = 1/2, 1, 5/2. Symbols are quantum MC data ^ 
for 5* = 1/2 (open circles), experimental data (filled symbols) 
for several real compounds (see FigJlClfor details), and the 
classical MC data (upwards triangles |E0|, downwards trian- 
gles jlS^, circles from our new simulations) used to construct 
rf 1 (fi), shown as dotted line; the open diamond is the exact 
result for the Ising model; see text for the dashed line. 



obtained by using the PQSCHA with the modified Low 
Coupling Approximation, as defined in Ref. Q , whose re- 
sults are shown by the dashed line. This modified version 
of the method has the disadvantage, however, of not being 
capable to cover the whole (—1,1) /^-interval, because of 
the loss of solution of its self-consistent equations when 
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Fig. 10. Ising critical temperature for the antiferromagnetic 
easy-axis model, near the isotropic limit. Curves and open 
symbols are as in Figjs^; Filled symbols are experimental 
data for the compounds YBa2Cu306.1 (downwards triangle), 
K2NiF4 (cross), Rb2Ni F4 (upwards triangle), Rh^MnCU (cir- 
cle), Rb2MnF4 (diamond); see text for references. 



moving towards the isotropic limit (reason for the dashed 
curve in Fig. ^|to be interrupted). 

On the other hand, the almost isotropic region, which 
is the most interesting from the experimental point of 
view, is well described by our results, even for S — 1/2, 
as seen in Fig. [To]. 



6 Conclusions 

The results found by the PQSCHA for the 2d Heisenberg 
antiferromagnets with easy-plane or easy-axis anisotropy, 
show that quantum fluctuations introduce temperature, 
anisotropy, and spin dependent renormalizations of both 
the energy scale and the anisotropy of the model, conse- 
quently reducing the temperature at which a phase tran- 
sition, of the BKT or Ising type, is expected. Such reduc- 
tion, however, is never seen to push the critical temper- 
ature down to zero, so that the whole class of quantum 
models described by Eq. (Q) with ^ |A| turns out to be 
characterized by a finite temperature transition. 

In Fig|ll] a global picture of the situation is visualized 
by the t c curves for different spin values, versus the ra- 
tio A//Lt: the central region corresponds to the easy-plane 
case, while the two external branches represents the easy- 
axis ferro- and antiferromagnetic cases. The t c reduction 
is seen to be more marked for small spin and in the ferro- 
magnetic case, thus leading to the conclusion that, should 
t c vanish for some critical v c (S), we expect | (5") | to be 
an increasing function of S which is larger in the antifer- 
romagnetic case, and such that its difference from unity is 
less than 10~ 2 . In the S = 1/2 case, despite being the one 
for which the PQSCHA is less precise, our picture corn- 



Fig. 11. Critical temperature vs X/fi for the class of models 
described by Eq. 1, with S = 1/2, 1, 5/2, and 00; symbols and 
curves as in Figs. 8 and ^[ 



pletely agrees with that proposed by Ding in Ref. [£9[ on 
the basis of quantum MC data. In the antiferromagnetic 
5 = 1/2 easy- axis case, on the other hand, Branco and 
De Sousa |3l[] predict, from a real space renormalization 
group analysis, the existence of a critical /i c (l/2) ~ —0.8 
at which t\ should vanish, and detect a reentrance in the 
\x dependence of such temperature; they found no such 
anomaly in the ferromagnetic sector. This scenario is in 
contrast not only with ours and Ding's results, but also 
with experimental data. 

Finally, from the experimental point of view, our re- 
sults show that anisotropies, even if as small as those 
measured in most real compounds, are fundamental in- 
gredients in determining the low-temperature behaviour 
of the materials; it is hence essential, in order to per- 
form a meaningful analysis of the experimental data, to 
have at least an estimate of the temperature at which a 
possible BKT or Ising-like phase transition should occur. 
Our approach allows a determination of such temperature 
for any given value of S and u, therefore representing an 
important tool in the experimental investigation of two- 
dimensional Heisenberg magnets. 
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